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Abstract 



HiggsBounds 2.0.0 is a computer code which tests both neutral and charged Higgs sectors 
. of arbitrary models against the current exclusion bounds from the Higgs searches at LEP 

and the Tevatron. As input, it requires a selection of model predictions, such as Higgs 
masses, branching ratios, effective couplings and total decay widths. HiggsBounds 2.0.0 
then uses the expected and observed topological cross section limits from the Higgs searches 
^ \ to determine whether a given parameter scenario of a model is excluded at the 95% C.L. 

by those searches. Version 2.0.0 represents a significant extension of the code since its first 
release (1.0.0). It includes now 28/53 LEP/Tevatron Higgs search analyses, compared to 
the 11/22 in the first release, of which many of the ones from the Tevatron are replaced 
by updates. As a major extension, the code allows now the predictions for (singly) charged 
Higgs bosons to be confronted with LEP and Tevatron searches. Furthermore, the newly 
included analyses contain LEP searches for neutral Higgs bosons {H) decaying invisibly or 
into (non flavour tagged) hadrons as well as decay-mode independent searches for neutral 
Higgs bosons, LEP searches via the production modes t'^t'H and hhH, and Tevatron 
searches via tiH. Also, all Tevatron results presented at the ICHEP'IO are included in 
version 2.0.0. As physics applications of HiggsBounds 2.0.0 we study the allowed Higgs 
mass range for model scenarios with invisible Higgs decays and we obtain exclusion results 
for the scalar sector of the Randall- Sundrum model using up-to-date LEP and Tevatron 
direct search results. 



1 Introduction 



The search for Higgs bosons is a major cornerstone of the physics programmes of 
past, present and future high energy colliders. The LEP and Tevatron experiments, 
in particular, have been able to turn the non-observation of Higgs bosons into con- 
straints on the Higgs sector, which can be very useful in restricting the available 
parameter space of particle physics models. At the LHC, discoveries involving colour 
charged particles (e.g. squarks, gluinos) may occur already before a Higgs signal gets 
established. The first LHC Higgs search results from ATLAS and CMS are expected 
to become public during the next months. It is important to note that limits from 
the Higgs searches at the LHC, the Tevatron and LEP will play a crucial role also 
for testing the validity of model interpretations of possible Higgs-like signals or other 
signals of new physics. 

The Higgs searches performed by the experimental collaborations usually take one of 
two forms. Dedicated analyses have been carried out in order to constrain some of 
the most popular models, such as the Standard Model (SM) [1] and various bench- 
mark scenarios in the Minimal Supersymmetric Standard Model (MSSM) [2]. In addi- 
tion, (essentially) model-independent limits on the cross sections of individual signal 
topologies (such as e~^e~ -> hiZ — bbZ) have been published. The former type of 
analyses include detailed knowledge of the overlap between the individual experimen- 
tal searches, and therefore have a high sensitivity, whereas the latter can be used to 
test a wide class of models. 

There are certain issues involved with the application of these experimental con- 
straints. The data is distributed over many different publications and the limits are 
given with a variety of normalisations. In the case of the Tevatron, the results are 
also frequently updated. Furthermore, care must be taken when using more than one 
experimental analysis to ensure the correct statistical interpretation of the result in 
terms of a confidence level (C.L.). One needs to take into account the underlying as- 
sumptions of the various experimental analyses in order to determine which analysis 
is applicable to which Higgs boson of a particular model. For example, if an analysis 
has been carried out assuming the validity of the SM, it will not be applicable to Higgs 
bosons whose various production and decay mechanisms do not show a proportionate 
similarity to the SM. 

The Fortran code HiggsBounds has been designed to facilitate the task of compar- 
ing Higgs sector predictions with existing exclusion limits, thus allowing the user to 
quickly and conveniently check a wide variety of models against the state-of-the-art 
results from Higgs searches. HiggsBounds 2.0.0 represents a significant extension 
of the code since its first release in February 2009 [3] and its latest release in De- 
cember 2009 (version 1.2.0). As a major extension, it allows now the predictions of 
(singly) charged Higgs sectors to be compared against current experimental limits 
from the LEP and Tevatron experiments. We have also significantly extended the 
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variety of neutral Higgs searches included. From the LEP experiments, we newly in- 
cluded searches where the Higgs boson is assumed to decay invisibly or into 77 or (non 
flavour tagged) hadrons. Also newly included are the decay-mode independent search 
of OPAL [4] and several searches involving the production process e~^e~ — ?■ ffH 
with / G {&, t}. From the Tevatron experiments, we newly included the analyses 
released publicly at the time of the ICHEP 2010 Conference, thus replacing many 
implemented analyses by their updated versions. Furthermore, we newly include the 
CDF and DO combination of their searches for the individual signal topology exclusive 
neutral Higgs boson production with decay into W~^W~ [5], which is applicable to a 
wide class of models beyond the SM, a search via the tiH production process of [6] 
and a search assuming the decay H — )■ 7Z [7] . A full list of the experimental analyses 
and the respective references included in HiggsBounds 2.0.0 is given in Sect. 4. 

In this release, we do not yet include charged Higgs search channels which involve 
the decay of charged Higgs bosons to neutral ones plus a W boson Hf — )■ h^W^. 
Neither are search results for doubly charged Higgs bosons included. Therefore, we 
will refer to the singly charged Higgs bosons below just as charged Higgs bosons. 

The rest of the paper is organised as follows. In Section 2, we give a brief description 
of the general approach of the program HiggsBounds2.0.0. Section 3 describes the 
model predictions required by the program as input in order to enable the application 
of Higgs search results. The three input options of HiggsBounds are described in 
detail, illustrating the input required from the user in each case and indicating how 
this input is further processed. Section 4 lists all analyses and their corresponding 
references implemented in version 2.0.0 of the program. Furthermore, the conditions 
implemented in the code which are responsible for determining the applicability of 
some analyses to certain Higgs bosons are described here. The complete operating 
instructions of HiggsBounds2.0.0 for the library of subroutines, the command-line 
and the online version with examples are given in Section 5. In Section 6, we present 
two new applications of the program. One application considers a SM-like Higgs boson 
where an extra Higgs decay channel to invisible particles is open; the other studies 
the scalar sector of the Randall- Sundrum model and shows, to our knowledge, for the 
first time the impact of Tevatron Higgs search results on this model. 



2 General Approach 

The following is a short, but self contained, description of the basic work-flow of 
HiggsBounds. A more detailed treatment, e.g. of the limitations of the methods 
employed here, can be found in the manual of HiggsBounds 1.0.0 [3]. 

The user provides the Higgs sector predictions of the model under consideration. For 
each neutral or charged Higgs boson Hi {i = 1, . . . ,nHO + nH±) in the model, this 
includes information on the Higgs mass, Higgs total decay width, Higgs and top quark 
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branching ratios and Higgs production cross sections, 



, rtot(-f^j) , BRmodcl(-f^j — ^ ■■■) , BRmodel(^ — ^ ■■■) , ™° „x , (1) 

and information whether each neutral Higgs boson is CP-odd, CP-even or is not a 
CP eigenstate. P stands for the Higgs production process. A full list of the required 
Higgs sector predictions is given in Sect. 3.1. Where it exists, a^^{P) is used as the 
reference production cross section. Variations on the input format (1) are offered, as 
described in detail in Sect. 3, which allow the user to specify a simpler set of input, 
provided certain basic approximations are valid. 

Note that the user may, to some extent, use the program also with incomplete model 
information by setting not provided cross section ratios, branching ratios or effective 
couplings to zero. However, the program may then consider some analyses not appli- 
cable, which would be taken into account if more complete model information had 
been provided. 

The HiggsBounds package includes sample files which demonstrate how the code can 
be used in conjunction with the widely used programs for predictions in the MSSM 
Higgs sector FeynHiggs [8-11] and CPsuperH [12,13]. The possibility to link Higgs- 
Bounds with the codes DarkSUSY [14], Fittino [15], MasterCode [16,17], 2HDMC [18] 
has already been demonstrated in various applications. 

A list of the experimental analyses currently included in HiggsBounds 2.0.0 can be 
generated by running the code. These include results from both LEP and the Teva- 
tron and consist of tables of expected cross section limits (based on Monte Carlo 
simulations with no signal) and observed ones at 95% C.L., with a variety of nor- 
malisations. The list mainly consists of analyses for which model-independent limits 
were published. However, we also include some dedicated analyses carried out for the 
case of the SM, or for Higgs bosons with certain CP properties. These analyses can 
only be applied to Higgs bosons which show the corresponding characteristics to a 
sufficiently high level of accuracy. 

We call the application of a Higgs search analysis to a particular Higgs boson (or two 
Higgs bosons, if the analysis involves two of them) of the model under study with 
particular mass(es) an "analysis application" X. ^ To each analysis application X 
corresponds a signal cross section prediction a{X) for the particular Higgs boson on 
which an upper limit is put by the analysis. For each X, HiggsBounds uses the input 
to calculate the relevant quantity Qmodei(-^) in terms of which the limit is given (i.e., 
a conveniently normalised cross section cr{X) times branching ratio). 

^ As an example, suppose a model with three neutral Higgs bosons {hi, /i2, /13) should be 
checked against two neutral Higgs search analyses Ai and A2, then there are six possible 
analysis applications X {Ai applied to hi, /12, /is and A2 applied to hi, /12, h^) that the 
model can be tested against. 
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The normalisation is carried out using SM predictions for Higgs boson production 
cross sections and decay branching ratios from HDECAY [19]version 3.4, the TEV4LHC 
Higgs Working Group [20] (for a full list of references, see Tab. 11), VFBSNLO [21], 
HJET [22-24] version 1.3 and dedicated calculations of our own [3]. 

The rationale behind the choice of SM normalisation is that virtually all Tevatron 
analyses implemented in this program use these predictions when normalising their 
cross section limits to SM quantities. Thus, describing deviations from the SM of a 
new model by using the SM normalisation the experimental analyses have chosen 
allows for the most accurate interpretation of the limits. By using SM predictions for 
Higgs production cross sections which deviate from the ones used in the experimen- 
tal analyses by a certain amount (be it because of a different loop order, different 
numerical values of input parameters, renormalisation scheme or choice of parton 
distribution functions), a deviation of the same relative size will be caused in the 
quantities Qmodci- The user should bear this in mind when interpreting the output of 
HiggsBounds. 

In order to ensure the correct statistical interpretation of the results, it is crucial 
to only consider the experimentally observed limit for one particular analysis ap- 
plication. Therefore, HiggsBounds must first determine Xq, which is defined as the 
analysis application with the highest statistical sensitivity for the model point under 
consideration. In order to do this, the program uses the tables of expected exper- 
imental limits to obtain a quantity Qcxpcc corresponding to each X. The analysis 
application with the largest value of Qmodci/Qexpec is chosen as Xq. 

HiggsBounds then determines a value for Qobs for this process Xq, using the appro- 
priate table of experimentally observed limits. If 



HiggsBounds returns as a result that this particular parameter point is excluded at 



In order to use HiggsBounds, the narrow- width approximation must be valid for 
each Higgs boson described in the input, because the experimental exclusion bounds 
currently utilised in the program have all been obtained under this assumption. Hig- 
gsBounds can be used with models which do not change the signature of the back- 
ground processes considerably ^ . In addition, models should not significantly change 

^ Note that if we had compared predicted cross sections directly against the experimentally 
measured limits for all available search channels and had excluded the model if at least one 
of them excludes it at 95 % C.L., the derived constraint would in general not correspond 
to a constraint at 95 % C.L. The maximally 5 % probability of each individual comparison 
of Qmodei and Qohs to yield a false exclusion would yield an overall probability for false 
exclusion which is higher than 5 %. 

^ However, new physics models which show strong deviations from the SM in the back- 




(2) 



95 % C.L. 2 
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the kinematical distributions of the signal cross section associated with X (e.g. r], 
Pt distributions of the final state particles) from that assumed in the corresponding 
analysis (for a full discussion see [3]). 

In addition, the LEP Higgs vector boson fusion (VBF) cross section should not be 
significantly enhanced compared to the LEP Higgsstrahlung cross section. (Models 
with a common ghiVV / 9hvv coupling for V = Z,W and no new mediating particles, 
such as the MSSM, satisfy this condition automatically.) 



3 Theoretical predictions required as input 

The user can choose between three different input formats, labelled by the variable 
whichinput. We describe here in detail what is required from the user for each of the 
three settings (hadr, part and ef fC). A summary can be found in Tab. 3, Sect. 5.1. 

3.1 The most general input form: whichinput=hadr 

The option whichinput=hadr requires the model input in the most general form. It 
involves specifying (at most) 

(1) masses for the neutral Higgs bosons {k = IjU^o) and singly charged Higgs 
bosons {j = l,nj^±), 

(2) Higgs total decay widths, 

rtot(^fc) 5 ^tot{H^), 

(3) whether the neutral Higgs boson is CP-even, CP-odd or mixed CP, 

(4) neutral Higgs branching ratios with SM equivalents (OP=ordinary particles i.e. 
particles which exist in the SM), 

BR„,odei(/ife -> OP) with OP = ss, cc, bb, r+r-, W^W-, ZZ, Z-f, -f-f, gg, 

(5) neutral Higgs branching ratios without SM equivalents 

BRmodoi(/ifc hihi),BR^odc\{hk invisible). 



ground processes to Higgs production are not usually considered in the literature, since this 
would often put them in conflict with SM electroweak precision data [25,26]. 
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(6) charged Higgs branching ratios to SM particles 

BR{H+ OP) with OP = cs, cb, r+z/,, 

(7) top quark branching ratios 

BR(t ^ W+b), BR(t ^ H+b), 

(8) normahsed cross section Ra{P) for LEP Higgs production process P 

e~^e~ hjZ, 
e~^e~ bbhj, 
e~^e~ —7- T^T~hj, 
e~^e~ — >■ hjhi, 

(9) normahsed cross section Ra{P) for Tevatron hadronic Higgs production process 
P 

PP hj, 

pp —J- bhj, 

PP — hjW, 

pp hjZ, 

pp — hj via VBF, 

pp — > tihj. 

Note that only a subset of these inputs may be required. For example, if only the 
LEP limits are required by the user, no Tevatron cross sections will need to be given 
as input (see Sect. 5). If the user only wishes to test the neutral Higgs sector, no 
input involving the charged Higgs sector or top decays will be required. 

The normalised cross section of a Higgs production process P is defined by 



T-i I n\ '^modcl 

O^vcdP) 

Where the SM equivalent exists, the reference cross section a^^fiP) of a process P 
involving a Higgs boson hk is crref(-P) = o-sm(-P)|, 



The only neutral Higgs production process without a SM equivalent considered in 
HiggsBounds up to now is P = e'^e~ — )■ hjhi. In this case, as in HiggsBounds 1.0.0 
[2], we choose a reference cross section for a fictitious production process of two scalar 
particles {H', H) with masses niH' = fUh. and mn = via a virtual Z exchange 
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with a standardised squared coupling constant: 



o2 



9wHz) — All ' 

where e denotes the electromagnetic coupling constant, and and the sine and 
cosine of the weak mixing angle, respectively. 

The cross section of this process in leading order is then completely determined by 
the Higgs masses and SM input and related to the SM Higgsstrahlung cross section 
via a simple phase space factor: 

are{{H' H production) = A (m///, mn, s) a^{mH) , (4) 

\{mH',mH,s) - ^'^^'^ 



^ _ {rria + nib 



2 
S 

21 



nia - nib) 



This reference cross section coincides with the MSSM tree level cross section for the 
process e~^e~ — )■ h^A^ if the Higgs mixing-angle dependent factor cos(/3 — a) is divided 
out of the tree level coupling and m^o and niho are chosen as niH' and nin- Therefore, 
i?o-(e^e~ — )■ h^A^) is simply given by cos^(/9 — a) in the MSSM (with real parameters) 
at tree level. 

For the process P = e^e^ — H^^HJ, the reference cross section is the cross section of 
the process e^e~ — > H^H^ in a 2 Higgs doublet model (e.g. the MSSM) at tree-level 
(i.e. s-channel 7 and Z exchange). This reference cross section depends solely on the 
mass of the charged Higgs boson and SM quantities. As a consequence, in the MSSM, 
R„{e+e- HfHj) = 1. 

The invisible branching ratio BRmodci(^fc invisible) is defined as the branching 
ratio of a neutral Higgs boson into particles which only appear in the detector as 
missing energy. Examples of invisible particles include neutralinos in the MSSM or 
majorons in supersymmetric models with spontaneous breaking of R-parity. 

The top quark branching ratios BR(t — )■ H'^'h) and BR(t — W^h) are also required 
as input, when the model is compared with the charged Higgs searches from the 
Tevatron. BR(t — >■ W^h) is required for analyses in order to determine when analyses 
which assume BR(t — )■ H'^h) -|-BR(t — )■ W^h) ~ 1 can be applied to the model point. 

In order to make it more convenient for users to calculate this input, the HiggsBounds 
library provides a series of Fortran functions which allow the user to access the predic- 
tions of certain SM quantities, including the SM Tevatron hadronic Higgs production 
cross sections, total decay width and branching ratios as a function of Higgs mass 
and r^^(t — )■ W^b) as function of the top mass. These are listed in Tab. 11 together 
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with the corresponding references. 



3.2 Specifying some cross sections at the partonic level: whichinput=part 



The input required for the option whichinput=part may be more convenient for the 
user to calculate than for the option whichinput=hadr, since most of the Tevatron 
cross section ratios can be specified at the partonic level rather than the hadronic 
level. 

The input option whichinput=part requires, at most, 

^l)-(8) as for whichinput=hadr 

(9) normalised cross section Ra{P) for Tevatron hadronic neutral Higgs production 
process P 

pp — hj via VBF 
pp —J- tihj 

(10) normalised cross section Rnm^ for Tevatron partonic neutral Higgs production 
process nm — hj + y, where y indicates a particle produced in association with 
the Higgs boson (or none if y is omitted): 

Rnin^ withnm =gg,bb, 
^nrn^^ , with um = ud, OS, 



Rnm^~ , with nm = du, sc, 
Rnin^ ^ with nm = dd, uu, ss, cc, bb, 
^h,+b,h,+-b^ with nm = bg, bg 



The normalised cross section Rnin^ for a Tevatron partonic neutral Higgs production 
process is defined by 



"model 

^ nm—^H+y 

calculated for a parton-system centre-of-mass energy squared s = so, where sq denotes 
the partonic production threshold sq = {mn + f^^y)^ (with my = in the case of single 
Higgs boson production). This requires that the dependence on s is (at most) mild. 

Internally, HiggsBounds uses the relation 

D ^p^_ V i^H+y ^SM{pP ^ nm ^ H + y) 
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to calculate the hadronic cross section ratios from the partonic cross section ratios. 
The ratios '^^^^^_^J^p^^^^^^ are provided within HiggsBounds. (For further explana- 
tion of the applicability of this approximation and details of how '^^CTsIf(ppTj^+i/)"^'' ^'^^ 
calculated, see [3]). 

The partonic cross section ratios R^m^ are much easier to calculate than hadronic 
cross section ratios in a wide range of models. In addition, it is often possible to make 
the approximation 

T3hj+W+ ^ T3hj+W- ^ T3hj+Z ^. T3hj + V 

for all nm, thus substantially reducing the number of partonic cross section ratios 
which need to be provided by the user from twelve to four. 

For instance in the MSSM with real parameters, this ratio R^^^^ can be calculated 
approximately from the normalised squared effective Higgs coupling to two Z bosons, 
i.e. idh'^zz / dnzzY (^^^ definition of this coupling, see Sect. 3.3). Similarly, RJ^ and 

R^g'lg'^^^'' can be taken to be approximately equal in the MSSM (and can also be 
calculated approximately from effective couplings: see Eq. (23)). 



3.3 The effective coupling approximation: whichinput=ef f C 



The option whichinput=ef f C allows the user to specify a greatly reduced number of 
input parameters. It involves specifying (at most) 



(1) masses for the neutral Higgs bosons h^. [k = l,nHo) and singly charged Higgs 
bosons (j = l,njj±), 

j 

(2) Higgs total decay widths, 

J- tot 

(3) normalised squared scalar and pseudoscalar effective couplings to fermions 



„modcl \ ^ / model ^ ^ 
^^.fefc(OP) I / ^P.fefc(OP) 

9h{op) J \ dnfop) 



OP = ss, CO, bb, tt, T~^T 



(4) normalised squared effective couplings to bosons 

' „model \ ^ / „model \ ^ 

A-W^] , OP = W-W~,ZZ,Z,,,,,gg 
9h'hz I \ 9h{ov) I 
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(5) neutral Higgs branching ratios without SM equivalents, charged Higgs branching 
ratios to SM particles and top quark branching ratios as before 

From this input, the input corresponding to the option part is calculated. 

We define the scalar and pseudoscalar Higgs coupling to fermions in the usual way, 
via the Feynman rule for the couphng of a generic neutral Higgs boson h to fermions: 

G{hff) = z{g,l + ig,^,) , (8) 

where Qs and Qp are real-valued scalar and pseudoscalar coupling constants respec- 
tively, and 1 and 75 are the usual matrices in Dirac space. A scalar particle, like the 
SM Higgs boson, has gp = Q and a pseudoscalar particle has (7^ = 0. 

Where it exists, the reference couplings squared (fifffYop))^ are the SM tree- level equiv- 
alents: 



= [^'^)\ (9) 
{9''^wwf=i-mw)\ (10) 




where mw and rrif denote the masses of the W boson and fermion /, respectively. 
idwHzY is defined in Eq. (3) above. 

The effective couplings (^7™°(OT)/fi'lf(OP))^' "^i^h OP = Zj, 77, are loop-induced. They 
can be defined via 



' „ model \ ^ 

^fefc(OP) I ^ 



(12) 



There is a choice of definition of the Higgs-gluon-gluon effective coupling squared 
i9h'°(pp)/9Hgg)^- '^^^ either be defined via decay widths: 



( model \ ^ pmodeU™, \ 

yHgg / H^OP\"''H )\mH=mh^ 

or via partonic cross sections: 



' „model \ 2 



„SM / "99 
9Hgg 



K^- (14) 
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In general, the input option whichinput=ef f C should only be used when both defi- 

-, model 



nitions result in very similar values for (fi'/j^(op) 



IdliggY- However, in certain circum- 



stances, this condition can be relaxed. For example, in a model in which the LEP 
searches for Higgs bosons decaying into hadrons are not relevant, the normalised 
gluon-gluon-Higgs effective coupling can be defined solely by Eq. (14). Conversely, if 
the Tevatron gluon fusion Higgs production mechanism is not relevant, the normalised 
gluon-gluon-Higgs effective coupling can be defined solely by Eq. (13). 



The normalised LEP cross sections are calculated from the normalised effective cou- 
plings using the relations: 



i?,(e+e- ^ hkZ) 



RJe+e- ^ bbh^^ 



R^{e+e- ^ r+r-/i^^^"^") 



/ „model \ ^ 

/ QhkZ" 




/ „modcl \ 2 
CPevenN _ / ^s.h^^T+T- 




(15) 
(16) 

(17) 



(19) 
(20) 



The analyses currently in HiggsBounds which use o"modei(e^e — )■ 66/ifc) or crmodei(e"'"e — )■ 
T^T~hk) apply only to Higgs bosons which are CP eigenstates. 

The partonic Tevatron cross section ratios are calculated by 



J^k j^hkb,hi^b 

bb bg,bg 



qq' ^q'q 



^q"q" 



model \ 2 
^h^gg 




(„model \ 2 
9hkWW \ 
9hww J 

C„ model \ ^ 
3hkZZ \ 
gfTzz) ' 



(21) 

(22) 

(23) 
(24) 



where (g, q') G {{u, d), (c, s)} and q" G {u, d, c, s, b}. 
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The tt CP-even Higgs hadronic Tevatron cross section ratio is obtained using 



Ra{pp tth)^ 



CP even \ 



^ „modej \ ■ 
9s,hktt \ 



(25) 



The normahsed hadronic cross sections for hk production via VBF is calculated using 
the approximate relation: 



(„model \ 2 / model \ ' 

9hww / V Ohzz / 



where the numbers 



R 



WW 
VBF 



R 



ZZ 
VBF 



c"sm(pp — ^ H via WW fusion) 

osuipp H via VBF) 
c"sm(pp H via ZZ fusion) 
csuipp — > H via VBF) 



77% 



23%. 



(26) 

(27) 
(28) 



have been calculated for pp colhsions with 1.96 TeV centre-of-mass energy using 
VBFNLO [21]. Including the Higgs mass dependence in the range 70 GeV < mn < 
300 GeV would change these proportions by less than 1%. Interference effects also 
affect the result by less than 1% (cf. [27]). (See Ref. [3] for a more detailed discussion.). 

Within HiggsBounds, the neutral Higgs decay widths to ordinary particles in the 
effective coupling input option are then calculated by 



pmodcl / 
hk-^OPV"'hk 



^modcl 



9hk{OP) \ T^SM 



SM 

y//(op) 



^H-^Opi^n] 



mH=mh, 



for OP = W+W-, ZZ, Z-f, -f-f, gg, 



^modcl 



9s,hk(OP) 

SM 

yH{op) 



-.model 



+ 



SM 
y//(OP) 



for OP = ss, cc, 66, t'^t" 



■pSM 



H^opy^H, 



(29) 
(30) 



with f3f{mh^] 



Am'j / . These are converted to branching ratios via 



T^pmodel / ^ _ ^h^^Opj^hk) 

aKy^^^Qp[mhJ - — - — 77-^^ — (31) 

Note that this means that it is especially important for the user to give an accurate 
rtot(^fc) when using this input option. 

Using effective couplings as input, the CP properties of the neutral Higgs bosons can 
be inferred, at least as far as the signal properties considered here are concerned. In 



13 



the present implementation, we set the CP value of a Higgs in the following way: 



For the implemented set of analyses, the limits on topological cross sections are 
applicable with high accuracy to models where the Higgs sector can be faithfully 
parametrised with effective couplings (See Ref. [3] for details). 

In summary, users may find it convenient to use the input option whichinput=ef f C 
for cases in which Eq. (6), Eqs. (15) - (26), and Eqs. (29) and (30) are valid to a 
satisfactory level of accuracy. Of course, if not all of the analyses are relevant to the 
model, some of these relations do not need to hold. We have already discussed that, 
in the case of a model in which the LEP searches for Higgs bosons decaying into 
hadrons are not relevant, Eq. (13) does not need to be valid. 



4 Analyses in HiggsBounds 2.0.0 

HiggsBounds contains many of the available observed and expected exclusion limits at 
95% CL. (preliminary and final results) from LEP and the Tevatron as data tables 
which are read in by the program during start-up. In order to provide values for 
Qobs(-^) and Qexpcci^) for continuous Higgs mass values, the program interpolates 
Q- values linearly between neighbouring Higgs mass points. We aim to include not only 
the most recent analysis in each channel by each collaboration (or analysis combining 
results from more than one collaboration) but also the most recent result appearing 
in a paper submitted to the arXiv. For simplicity we use the label 'published' for 
those results (since results submitted to the arXiv are usually being prepared for 
publication in a journal or conference proceedings). 

In HiggsBounds 2.0.0, in total 81 Higgs search analyses have been implemented 
consisting of 28 from LEP and 53 from the Tevatron. While many types of analyses 
have been added for the first time, several Tevatron analyses included in HiggsBounds 
1.2.0 have been replaced by updated ones based on more data. For completeness and 
as a reference for users, a full list of the implemented analyses follows. 






CP(/ifc) -> 



if neither of the above conditions is met. 
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• LEP neutral Higgs analyses (considering the final states): 



hkZ, hk bb or tt [2], h^Z, ht hihi, hi bb [2], 

hkZ, hk anything [4], h^Z, hj, -> hihi,hi^ tt [2], 

hkZ, hk — invisible [28-31], h^hi, hk,i bb [2], 




hkZ, hk hadrons"^ , hkhi, hk — )■ hihi, hi — )■ bb [2], 

bbhk — )■ bbbb, hk CP even or odd, [33], hkhi, hk — )■ hihi, hi — )■ tt [2], 



bbhk — )■ bbTT, hk CP even or odd, [33,34], hkhi, hk — )■ bb, hi — )■ tt [2], 
rr/ifc — rrrr, hk CP even or odd, [33]. 

• Tevatron single topology neutral Higgs analyses (considering the final states): 



bhk bTT [56-58]. 

• Tevatron SM Higgs combined analyses considering the final states (schematic): 
Vhk ^bb + Ef"-^ with V E {W, Z} [59-62], 

Vhk -> VVV ^ l^l^ + X with / e {e,/x} and V G {W, Z} [63], 

hk + X ^ WW + X [64-70], 

hk + X ^TT + X [71-73], 

hk + X ^bb + X [74], 

hk + X^^^ + X [75-77], 

/ife + X and various Higgs decays [78-85]. 

• LEP/Tevatron charged Higgs analyses: 
e+e- ^ H+H^ 4 jets [86,87], 

e+e" H+H- tutu [87], 

pp ^ti,t^ H+b, H+ cs (& c.c.) [88,89], 

pp ^tt,t-> H+b, H+ Tu (& c.c.) [88]. 

^ Combination based on [35-38] 



Zhk llbb [39-41] , 
Whk -> lubb [43-46], 
Whk 3W, [51,52], 
bhk 36 jets [53-55] 



single hk WW [42,5] , 
single hk tt [47-50], 
single hk-^ Z'j, [7], 
tihk tibb [6], 
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Internally, HiggsBounds uses a number of predictions for SM quantities for the Higgs 
sector[19-24, 90-107] to convert between experimental limits with different normali- 
sations. 



Some of the analyses have been performed under certain model assumptions. Typical 
assumptions are, likeness to the SM, or models which fulfil BR(t — W~^b) + BR(t — )■ 
iJj^fo) = 1, or Higgs bosons with certain CP properties. In these cases, HiggsBounds 
determines for a given model scenario which analyses can be applied to which Higgs 
boson. 

In order to decide whether a model scenario is sufficiently 'SM-like' to be compared 
with an analysis carried out under SM assumptions, we use the following test. Firstly, 
for each of the A'^cs distinct production cross sections crmodci{Pi{h)) and the Nbr 
distinct decay branching ratios BR(/i — )■ F^) which appear in the list of M signal 
topologies in that particular analysis, we determine the normalised mean value and 
the deviation from the mean: 

^ = E ' ^ = TT- , (32) 

iVcs i=i ^Vbr 

5si = Si-s, 6bk = bk-b, (33) 

.,1 O"modcl(-Pj(^)) , BRmodcl(^ — ^ -^fc) /^.x 

"""" " = asuiP.m ' = BRsM(ff F.) ■ <^^' 



The parameter point is considered 'SM-like' for this particular analysis if the relative 
deviations of all individual production and decay combinations (z, k) from the mean 
sb stay below a preset bound e: 



max 

i,k 



Ssi , 5bk , SsiSbk 



s b sb 



<e. (35) 



By default, e = 2%, i.e. the predictions for the different topological cross sections, 
normalised to the SM values, are required to result in the same common scale factor 
with at least 2% accuracy. This method is conservative and errs on the side of caution, 
but can occasionally be overly restrictive. 

After ensuring that a certain parameter point passes the SM-likeness test for this 
analysis, the normalised theoretical cross section Qmodei is calculated simply by 



Qmodci = sb. (36) 

If two or more neutral Higgs bosons have similar masses, then HiggsBounds allows 
the possibility of adding up their cross sections. If the mass difference between two 
neutral Higgs bosons is less than delta_Mh_LEP (or delta_Mh_TEV), then the full LEP 
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(or Tevatron) cross sections will be added. The normalisation factor is calculated at 
the average mass of the Higgs bosons involved. Note that: 

delta_Mh_LEP and delta_Mh_TEV are set to zero by default. 

To switch this feature on, the user should change these variables in the code and 
recompile (see the readme files about where to find these variables in the code). 
This should only be done if the user is confident that interference effects are small. 
We would recommend that delta_Mh_LEP < 2 GeV and delta_Mh_TEV < 10 GeV. 
For analyses which require a SM-likeness test, the full cross sections are not added, 
regardless of the values of delta_Mh_LEP and delta_Mh_TEV. 

This feature can be very useful in, for example the rn^^^ scenario of the MSSM [108], 
where a CP even Higgs boson and the CP odd Higgs boson can have very similar 
masses but, naturally, they do not mix. In this case, the cross sections are commonly 
added, such as for example in [53] and [54]. 



5 New HiggsBounds Operating Instructions 

Features which are new in version 2.0.0 will be labelled with the symbol: *. 
There are three formats in which the program HiggsBounds can be used: 

• Library of subroutines 

• Command-line version 

• Online version 

The most widely applicable format of HiggsBounds is the command-line version, since 
this reads all the model data from text files and thus this model data can be generated 
using any package the user wishes. The library of HiggsBounds subroutines allows 
HiggsBounds to be called within other programs. If the user just wishes to check a 
few parameter points, the online version provides quick access to all the functionality 
of HiggsBounds, without the need to install the code. 

The HiggsBounds code, the online version and documentation can all be found at 
the URL http://projects.hepforge.org/higgsbounds/ . ^ 

The HiggsBounds code is provided in either Fortran 77 or Fortran 90. Both codes 
provide exactly the same functionality and have exactly the same operating instruc- 
tions. In fact, the Fortran 77 and Fortran 90 versions of the HiggsBounds subroutines 
can even be called within codes written in Fortran 90 and Fortran 77, respectively. 
Both codes have also been tested with a variety of Fortran compilers, including the 

^ The former website www . ippp . dur . ac . uk/HiggsBounds should redirect to the new one. 
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free gnu compilers which accompany most Linux distributions. Therefore, the user 
may install either code and the difference will only be apparent if the user wishes to 
examine the structure of the code ^ . 

The library of subroutines, the command-line version and the online version share 
a common set of features, which we will describe first. We will then give operating 
instructions for each of these three HiggsBounds formats individually. 

5. 1 Common features: Input 

HiggsBounds requires five types of input: 

• the number of neutral Higgs bosons in the model under study (nHzero) 

• the number of singly, positively charged Higgs bosons in the model under study 
(nHplus) 

• the set of experimental analyses which should be considered (whichanalyses) 

• the theoretical predictions of the model under study (a set of input arrays) 

• the format of these theoretical predictions (whichinput) 

Tab. 1 contains further information on the variable nHzero and nHplus, and the 
possible values of whichanalyses are described in Tab. 2. Note that the option 
whichanalyses=' singH' should only be used if neither processes involving Higgs 
pair production at LEP nor processes involving the hj — >■ hihi decay are relevant. 
However, if these conditions are met, this option can save significantly on computing 
time. 

HiggsBounds expects the theoretical input to be in one of three formats, labelled by 
the variable whichinput. These formats were described in detail in Sect. 3 and are 
briefly summarised in Tab. 3. In Tab. 4, Tab. 5, Tab. 6 and Tab. 7 we assign names 
to all of the possible input arrays (each array is defined in terms of the notation used 
in Sect. 3). These names will prove useful when we describe the input requirements 
of each version of HiggsBounds individually. 

nHzero nHplus* (integer) 

0-9 0-9 This is the default range. 

It can easily be extended by the user if required. 

Table 1 

The possible values of the variable nHzero, which labels the number of neutral Higgs bosons 
in the model under study and the variable nHplus, which labels the number of singly, 
positively charged Higgs. 



" The maintenance of two separate codes is primarily intended to provide an efficient way 
for the authors to confirm that each update is free from implementation errors. 
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whichanalyses (character (LEN=5) ) 



LandT 


both LEP and Tevatron analyses 


onlyL 


only LEP analyses 


onlyT 


only Tevatron analyses 


onlyP* 


only published analyses (defined as analyses with an arXiv number 


singH 


only analyses for processes involving one Higgs boson 


Table 2 





The possible values of the variable whichanalyses, which indicates which subset of exper- 
imental analyses will be considered by HiggsBounds. 



whichinput (character (LEN=4)) 

ef f C Masses, total decay widths, 

ratios of effective couplings squared, some branching ratios. 

part Masses, total decay widths, ratios of LEP cross sections, 

mainly ratios of partonic Tevatron cross sections, branching ratios. 

hadr Masses, total decay widths, ratios of LEP cross sections, 

ratios of hadronic Tevatron cross sections, branching ratios. 

Table 3 

The possible values of the variable whichinput, which indicates the format of the theoretical 
predictions provided by the user for the neutral Higgs sector. (See Sect. 3 for a more detailed 
description of each of these settings). 
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input arrays 



(double precision) 



Mh(nHzero) 




in 


GeV 


MhGammaTot (nHzero) 


^ tot (hi) 


in 


GeV 


Mhplus(nHplus) * 




in 


GeV 


MhplusGammaTot (nHplus) * 


rtot(i/f) 


in 


GeV 



g2hjss_s (nHzero) * 

g2hjcc_s (nHzero) * 
g2hjbb_s (nHzero) * 
g2hjtoptop_s (nHzero) 
g2hjtautau_s (nHzero) 

g2hjss_p (nHzero) * 

g2hjcc_p (nHzero) * 
g2hjbb_p (nHzero) * 
g2hjtoptop_p (nHzero) 
g2hjtautau_p (nHzero) 

g2hj WW (nHzero) 

g2hjZZ (nHzero) 
g2hjZga (nHzero) * 
g2h j gaga (nHzero ) 
g2hjgg (nHzero) 



s,hj{OP) 

3M 

^//(OP) 



OP 



P,hj {OP) 

THM 

yH(OP) 



OP 



^model \ 2 
fhj{OP) 

THH 

i'H(OP) 



OP 



g2hj hiZ (nHzero , nHzero) 



^ model \ 2~ 



ss 

cc 
hi 
tt 

ss 

cc 
hi 
tt 

w+w- 
zz 

11 

99 



Table 4 

Input arrays for model predictions for effective normalised squared couplings recognised by 
HiggsBounds. The size of each array is given in brackets in the first column. See Sect. 3 for 
the description of the notation used in the second column. 
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input arrays cont. (integer) 



CP_value (nHzero) *, where CP_value(i) =-1 if ith neutral Higgs is CP-odd 

= if zth neutral Higgs has mixed CP 
= 1 if ith neutral Higgs is CP-even 

Table 5 

Input array to specify the CP properties of each of the neutral Higgs bosons predicted by 
the model. The size of the array is given in brackets in the first column. 



input arrays cont. 


(double precision) 




BR_hjss(nHzero) * 




OPl OP — 


ss 


BRJijcc(nHzero) * 






cc 


BRJijbb(nHzero) 






bb 


BR_hjtautau(nHzero) 








BRJijWW(nHzero) 






w+w- 


BRJijZZ(nHzero) 






zz 


BR_hjZga(nHzero) * 






Zl 


BRJi j gaga (nHzer o ) 






77 


BRJijgg(nHzero) * 






99 


BRJij invisible (nHzero) * 


BR{hj - 


-)• invisible) 




BRJijhihi (nHzero , nHzero) 


BR{hj - 


-> hihi) 




BR_tWpb * 


BR(t 


W+b) 




BR_tHpjb(nHplus) * 


BR(t 


Hp) 




BR_Hpjcs(nHplus) * 


BR{H+ 


OP), OP = 


cs 


BR_Hpjcb(nHplus) * 






cb 


BR_Hpjtaunu(nHplus) * 









Table 6 

Input arrays for model predictions for branching ratios recognised by HiggsBounds. The 
size of each array is given in brackets in the first column. See Sect. 3 for the description of 
the notation used in the second column. The elements of BRJijhihi are ordered such that 
BR_hjhihi(j ,i)= ^ hihi). 
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input arrays cont. 



(double precision) 



CS_lep_hj Z_ratio (nHzero) 


RAP)-, 


P = 




—7- hjZ 


CS_lep_bbhj _ratio (nHzero , nHzero) * 








— >■ bbhj 


CS_lep_tautauhj _rat io (nHzero , nHzero) * 






e~^e~ 


— >■ T^T^hn 
1 11/ J 


CS_lepJijhi_ratio (nHzero , nHzero) 








— > hnhi 
J ^ 


nS 1 pn Hni Hm "i rat i o CnH^pro") * 










CS_tev_ppJij _ratio (nHzero) 








hi 




V VJ LJ. 1 VJ L CL L> J. W \ J.J.11^ C J. W / 






TiT) — ^ 


bh ■ 


CS_tev_ppJij W_ratio (nHzero) 






Vf) —7- 


3 


CS_'tev_pp_h.j Z_ratio (nHzero) 






T)f) — y 




CS_tev_pp_vbf _ratio (nHzero) 








/ij viaVBF 


CS_tev_pp_tthj _ratio (nHzero) * 






vv — y 


tthj 


CS_tev_ggJij _ratio (nHzero) 




nm = 


QQ 




CS_tev_bbJij _ratio (nHzero) 






bb 




CS_tev_udJij Wp_ratio (nHzero) 




nm = 


ud 




CS_tev_cs_hj Wp_ratio (nHzero) 






cs 




nFi tp^T" nH Vi i Utti ral"! o rnH^prn') 

WkJ U C V LIV^ LJ. 1 VVIii L CL O J. W V, W / 


j^h,+W- 


JlVfl — 


dvj 




CS_tev_cs_h.j Wiii_ratio (nHzero) 






sc 




CS_tev_ddJij Z_ratio (nHzero) 


h,+Z 


nm = 


dd 




CS_tev_uuJij Z_ratio (nHzero) 






uu 




CS_tev_ss Jij Z_ratio (nHzero) 






ss 




CS_tev_ccJij Z_ratio (nHzero) 






cc 




CS_tev_bbJij Z_ratio (nHzero) 






bb 




CS_tev_bgJijb_ratio (nHzero) 


jjhj+b,hj+b 


nm = 


bg,bg 





Table 7 

Input arrays for model predictions for cross section ratios recognised by HiggsBounds. The 
size of each array is given in brackets in the first column. The LEP or hadronic Tevatron 
cross section ratios Ra{P) are defined Eq. (3), and the partonic Tevatron cross section 
ratios Rnm^ are defined in Eq. (5). 
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5.2 Common features: Output 



HiggsBounds provides the user with four types of output: 

• whether the parameter point is excluded at the 95% CL. or not (HBresult) 

• the reference number of the analysis application (Xq) with the highest statistical 
sensitivity (chan). 

• the number of Higgs bosons which have contributed to the theoretical rate for the 
corresponding process (ncombined) 

• the ratio of the theoretical rate Qmodei to the observed limit Qobs for this process 
(obsratio). 

Tab. 8 shows the possible values of HBresult and obsratio, which are complemen- 
tary. Tab. 9 and Tab. 10 gives information on chan and ncombined respectively. If 
the library of subroutines or the command-line versions are used, the key associating 
the reference numbers with the analysis applications is written in the file Key.dat. 
In the online version, this information appears on the screen. 



HBresult obsratio 

(integer) (double precision) 

> 1.0 parameter point is excluded 

1 < 1.0 parameter point is not excluded 

-1 < 0.0 invalid parameter set 

Table 8 

The possible values of the output variables HBresult and obsratio, which indicate whether 
a parameter point has been excluded at the 95% C.L. by the experimental results under 
consideration. 



chan (integer) 

0- [t^ of considered analyses] See the file Key . dat for the definition of each 

reference number. Key.dat is automatically 
generated when either the command line or the 
subroutine version of HiggsBounds are used. 

Table 9 

Further information about the output variable chan, which stores the reference number of 
the analysis application with the highest statistical sensitivity. (0= no process applies ) 
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ncombined (integer) 



1-nHzero Number of neutral (or singly, positively charged) Higgs bosons 
which have contributed to the theoretical rate for this process. 

Table 10 

Further information about the output variable ncombined. 
5. 3 Library of subroutines 

Installation 

The HiggsBounds code can be compiled to form a library of subroutines using the 
following commands: 

. / configure 
make libHB 

A program for which the HiggsBounds subroutines should be used can be compiled 
and linked to the library by adding -L<HBpath> -IHB to the command line, for 
example, 

gfortran myprog.fQO -o myprog -L<HBpath> -IHB 
where <HBpath> is the location of the HiggsBounds library. 

The HiggsBounds subroutines make use of the file handles 10, 11, 44, 45 and 87 
(i.e. users should avoid using these file handles in the program which calls these 
subroutines.). 

Subroutine initialize_HiggsBounds ^ 

The subroutine initialize_HiggsBounds must be called before any other Higgs- 
Bounds subroutine. It performs some preparatory operations such as reading in the 
tables of data. It is called as: 

call initialize_HiggsBounds (nHzero , nHplus, whichanalyses) 

When using the subroutines in another code, the subroutine initial ize_HiggsBounds 
must be called only once, before any other HiggsBounds subroutine is called. 

If the user does not wish to test the neutral Higgs sector with HiggsBounds, then 
he/she should set nHzero=0. Similarly, if the user does not wish to test the charged 
Higgs sector with HiggsBounds, then he/she should set nHplus=0. 
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Subroutines HiggsBounds_neutral_input_ef f C *, 

HiggsBounds_neutral_input_part *anc? HiggsBounds_neutral_input_hadr * 

The neutral Higgs sector input can be passed to HiggsBounds using one of the sub- 
routines HiggsBounds_neutral_input_ef f C, HiggsBounds_neutral_input_part or 
HiggsBounds_neutral_input_hadr. They set the value of whichinput to be ef fC, 
part and hadr respectively and therefore require different arguments. The use of one 
of these subroutines is only required if nHzero is non-zero (recall that nHzero is set 
in subroutine initialize_HiggsBounds). These subroutines are called as: 

call HiggsBounds_neutral_input_ef f C(Mh,MhGainmaTot , 

& g2h j ss_s , g2hj ss_p , g2hj cc_s , g2hj cc_p , g2hjbb_s , g2hj bb_p , 
& g2h j toptop_s , g2hj toptop_p , g2hj tautau_s , g2hj tautau_p , 
& g2h j WW , g2h j ZZ , g2h j Zga , g2h j gaga , g2h j gg , 
& g2hjhiZ,BR_hjinvisible,BR_hjhihi ) 

call HiggsBounds_neutral_input_part (MhjMhGammaTot , CP_value, 



& 


cs. 


.lep_hj Z_ratio , CS_lep_bbhj _ratio , CS_lep_tautauhj _ratio , 


& 


cs. 


.lep_hjhi_ratio , 


& 


cs. 


.tev_gg_hj _ratio , CS_tev_bb_hj_ratio , 


& 


cs. 


.tev_bg_hjb_ratio , 


& 


cs. 


.tev_ud_hjWp_ratio ,CS_tev_cs_hjWp_ratio , 


& 


cs. 


.tev_ud_hj Wm_ratio , CS_tev_cs_hj Wm.ratio , 


& 


cs. 


.tev_dd_hjZ_ratio, CS_tev_uu_hjZ_ratio , 


& 


cs. 


.tev_ss_hjZ_ratio, CS_tev_cc_hjZ_ratio , 


& 


cs. 


.tev_bb_hjZ_ratio, 


& 


cs. 


.tev_pp_vbf _ratio , CS_tev_pp_tthj_ratio , 


& 


BR. 


.h j s s , BR_h j c c , BR_h j bb , BR_h j t aut au , 


& 


BR. 


.h j WW , BR_h j ZZ , BR_h j Zga , BR_h j gaga , BR_h j gg , BR_h j invi s ible , 


& 


BR. 


.hjhihi ) 



call HiggsBounds_neutral_input_hadr (MhjMhGammaTot , CP_value, 

& CS_lep_hj Z_ratio , CS_lep_bbhj _ratio , CS_lep_tautauhj _ratio , 

& CS_lep_hjhi_ratio, 

& CS_tev_pp_hj_ratio , CS_tev_pp_hjb_ratio , 

& CS_tev_pp_hjW_ratio , CS_tev_pp_hjZ_ratio , 

& CS_tev_pp_vbf _ratio , CS_tev_pp_tthj_ratio , 

& BR_h j s s , BR_h j c c , BR_h j bb , BR_h j t aut au , 

& BR_h j WW , BR_h j ZZ , BR_h j Zga , BR_h j gaga , BR_h j gg , BR_h j invi s ible , 

& BR.hjhihi ) 

Each of these arguments must be supplied. However, if a branching ratio, effective 
coupling or cross section is believed to be irrelevant, the corresponding array may 
be filled with zeros. This will ensure that the value of Qmodei for processes involving 
this quantity will also be zero. For example, in the MSSM, scenarios exist where the 
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decay hj 77 will certainly not be the process with the highest statistical sensitivity 
of Qmodci/Qexpec and, Consequently, it may be convenient to set the arrays g2hjgaga 
and BR_hjgaga to zero for simplicity in such a case. Note that if a quantity used in 
a SM-likeness test is set to zero, the model point will fail that SM-likeness test, even 
if the corresponding SM quantity is very small. 

Also, depending on the value given for whichanalyses, some of the input arrays 
will be ignored within HiggsBounds. For example, if whichanalyses='onlyT', the 
branching ratio for the Higgs cascade decay hj — )■ hihi will not be relevant. Therefore, 
setting this array to zero will not affect the HiggsBounds results in this case. 



Subroutine HiggsBounds_charged_input * 

The subroutine HiggsBounds_charged_input gives the charged Higgs sector input to 
HiggsBounds. The use of this subroutine is only required if nHplus is non-zero (recall 
that nHplus is set in subroutine initialize_HiggsBounds). It is called as: 

call HiggsBounds_charged_input (MHplus , MHplusGammaTot 
& CS_lep_HpjHmi_ratio, 

& BR_tWpb,BR_tHpjb,BR_Hpjcs,BR_Hpjcb,BR_Hptaunu ) 



Subroutine run_HiggsBounds * 

This subroutine performs the main part of the HiggsBounds calculations. It is called 
as: 

call run_HiggsBounds ( HBresult , chan , 

& obsratio, ncombined ) 



Subroutine f inish_HiggsBounds 

The subroutine f inish_HiggsBounds should be called once at the end of the program, 
after all other HiggsBounds subroutines . It is called as: 

call f inish_HiggsBounds 



' In the Fortran 90 version of the code, the subroutine f inish_HiggsBounds is used to 
deallocate the allocatable arrays used within HiggsBounds. 
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Functions for SM branching ratios, total decay width and cross sections 

The HiggsBounds library also allows users access to the SM Higgs branching ra- 
tios, total decay width and production cross sections, which are used internally by 
HiggsBounds. We use SM Higgs branching ratios and total decay width from the 
program HDECAY 3.4 [19]. The SM hadronic cross sections have been obtained from 
the TEV4LHC Higgs Working Group [20] (see Tab. 11 for references to the original 
works) with the exception of the a^^{pp bg ^ bH) and a^^{pp -^bg^ bH) cross 
sections. The latter cross sections have been calculated with the program HJET 1.3 
[22-24] for a set of different cuts on the transverse momentum and pseudo-rapidity 
of the 6-quark, which are needed internally in order to apply correctly the results of 
some Tevatron analyses and which were not available from [20]. From this set, only 
the cross section without cuts is externally provided. 

Also included is a function for the SM top quark decay width into W^b in GeV, 
which depends on the top quark pole mass in GeV, given at next-to-leading-order, 
neglecting terms of order m^/m^, a] and (as/^r) M^/m^, as quoted in [109] (original 
reference: [110]). This is not used within HiggsBounds, but may prove useful when 
calculating BR(t H^b). 

Note that these functions are only valid within certain mass ranges, which are suffi- 
cient for the requirements of HiggsBounds. If a function is called with an argument 
outside this range, it returns a value of —1. 



Examples 

We have provided three example programs which demonstrate the use of the Higgs- 
Bounds subroutines. The first example relates to the Fourth Generation Model and is 
contained in the file example-SM_vs_4thGen.F. This program uses the HiggsBounds 
functions for the SM branching ratios and SM total decay width to calculate the 
Higgs decay width and the effective normalised squared couplings in the SM and 
a simple Fourth Generation Model. This information is then used as input for the 
subroutine HiggsBounds_neutral_input_ef f C, which is called once with SM input 
and once with Fourth Generation Model input. Once the HiggsBounds library has 
been compiled (using . /configure ; make libHB as described previously), the code 
example-SM_vs_4thGen.F can be compiled and run with the commands: 

gfortran example-SM_vs_4thGen.F -o example-SM_vs_4thGen \ 

-L<HBpath> -IHB 
. /example-SM_vs_4thGen 

where <HBpath> is the location of the HiggsBounds library. 
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function (double precision) 



C! -1 rum -i Vi/' MVi ^ 

oribcLninia._Li i^riiij 


J- tot ['^i) 


^^ Ql 


SMBRJiss (.Mnj 


rJK [H — > UFj, Ur= ss 


[19J 


SMBRJicc (.Mn; 


cc 


[19J 


SMBR_Hbb(,Mn; 


00 


[19J 


SMBRJltoptop (Mn; 


4-4- 

it 


ri nl 
[19J 




T 


[1 Ql 


briDHJiiiiuiiiu u^-ti^ 




Tt ol 


SMBR_HWW(.Mnj 




Fi nl 

[19J 


SMBR_HZZ(Mh) 


zz 


[19] 


SMBR_HZgam(Mh) 




[19] 


SMBR_Hgamgam(Mh) 


77 


[19] 


SMBR_Hgg(Mh) 


gg 


[19] 



SMGamma_tWpb(Mtop) * r^^{t ^W+b) [110] 



SMCS_tev_pp_qq_HW(Mh) cr^^(P), 
SMCS_tev_pp_qq_HZ (Mh) 
SMCS_tev_pp_gg_H (Mh) 
SMCS_tev_pp_bb_H (Mh) 
SMCS_tev_pp_vbf_H(Mh) 
SMCS_tev_pp_ttH(Mh) * 
SMCS_tev_pp_bg_Hb(Mh) a^^{pp^bg 
Table 11 

SM Higgs branching ratios, total decay widths in units of GeV and hadronic Tevatron cross 
sections in units of pb provided as functions by HiggsBounds, together with references. Each 
function takes a Higgs mass Mh (double precision) as its argument. Also included is the 
SM top quark decay width into W^b in GeV, which depends on the top quark pole mass 
in GeV. If a function is called with an argument outside its mass range, it returns a value 
of-1. 



P= pp^qq^ HW [98-100] 

pp^qq^ HZ [98-100] 

pp^gg^H [90-97] 

pp^bb^ H [101] 

pp-^ H via VBF [2 1 , 1 02-1 04] 

pp tin [105-107] 

Hb) + fjS^ (pp -^bg^ Hb) [22-24] 



The files HBwithFH.F and HBwithCPsuperH.f demonstrate the use of the subrou- 
tine version of HiggsBounds with the pubhcly available programs FeynHiggs [8-11] 
and CPsuperH [12,13], respectively. We refer the reader to the extensive comments 
contained within these example files for further details. 
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5-4 Command line version 



Installation 

In order to be able to call HiggsBounds from the command line, it should be compiled 
using the commands 

. / configure 
make 



Command line and input file format 

In the command-line usage of HiggsBounds, the arrays containing the theoretical 
model predictions are read from text files. The other options are specified in the 
command line, which is of the form: 

./HiggsBounds <whichanalyses> <whichinput> <nHzero> <nHplus> <prefix> 

The variable <pref ix> is a string which is added to the front of input and output 
file names and may include directory names or other identifying information. 

Tab. 12 and Tab. 13 describe the contents of each input file. Note that each input file 
should start with a line number. The input files should not contain any comments 
or blank lines. The line number identifies the predictions which belong to the same 
model parameter point in different files. The file BR_H_OP . dat * contains neutral Higgs 
branching ratios with a SM equivalent (i.e. ' OP '= 'ordinary particles'), whereas the 
file BR_H_NP . dat * contains neutral Higgs branching ratios without a SM equivalent 
(i.e. 'NP'='new particles') ^ . 

Care should be taken with the order of the array elements in the files. This is best 
illustrated by an example, where we will use uh = 3. The one dimensional arrays, 
e.g. Mh, should be given in the order 

Mh(l), Mh(2), Mh(3). 



In previous versions of HiggsBounds, this was equivalent to dividing the neutral Higgs 
branching ratios between a file containing one type of neutral Higgs and a file containing 
two types of neutral Higgs bosons. However, the input to HiggsBounds 2.0.0 requires the 
branching ratio of the Higgs to invisible particles, hence the new filenames. 
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file name data format 



MH_GaniniaTot . dat 


k, 


Mh, MhGammaTot 


MHplus_GammaTot . dat * 


k, 


Mhplus , MhplusGammaTot 


effC.dat W 


k, 


g2hj ss_s , g2hj ss_p , g2hj cc_s , g2hj cc_p 






g2h j bb_s , g2hj bb_p , g2hj toptop_s , g2h j toptop_p 






g2h j t aut au_s , g2h j t aut au_p , 






g2hjWW,g2hjZZ,g2hjZga, 






g2hjgaga,g2hjgg. 






some elements of g2hjhiZ 






(lower left triangle - see example) 


LEP_HZ_CS_ratios . dat 


k, 


CS_lepJij Z_ratio 


LEPJI_f f _CS jratios . dat * 


k, 


CS_lep_bbhj_ratio, CS_lep_tautauhj jratio 


LEP_2H_CS_ratios . dat 


k, 


some elements of CS_lepJijhi_ratio 






(lower left triangle - see example) 


LEP_HpHm_CS_ratios . dat * 


k, 


CS_lep_Hp j Hmj _rat io 


TEV_H_Ojet_partCS_ratios.dat 


k, 


CS_tev_ggJij _ratio , CS_tev_bbJij _ratio 


TEV_H_ljet_partCS_ratios.dat 


k CS_tev_bgJijb_ratio 


TEV_HW_partCS_ratios . dat 


k, 


CS_tev_udJijWp_ratio,CS_tev_csJijWp_ratio, 






CS_tev_udJij Wmjrat io,CS_tev_cs Jij Wm_rat io 


TEV_HZ_partCS_ratios . dat 


k, 


CS_tev_ddJijZ_ratio, CS_tev_uuJijZ_ratio , 






CS_tev_ssJijZ_ratio, CS_tev_ccJijZ_ratio , 






CS_tev_bbJij Z_ratio 


TEV_H_vbf JiadCSjratios . dat 


k, 


CS_tev_pp_vbf _ratio 


TEV_H_ttJiadCS_ratios . dat * 


k, 


CS_tev_pp_tthj jratio 


TEV_1H JiadCSjratios . dat W 


k, 


CS_tev_ppJij_ratio, CS_tev_pp_hjbjratio, 






CS_tev_ppJijW_ratio,CS_tev_ppJijZ_ratio , 






CS_tev_pp_vbf _ratio 






CS_tev_pp_tthj jratio 



Table 12 

Names and data format of all HiggsBounds input files (part I). The right column shows 
the order of the input data arrays within one line of the input file. For the order within 
the arrays, see example, k is the line number. Note that the arrays CS_tev_pp_vbf _ratio, 
CS_tev_pp_tthj_ratio each appear in two different input files. However, these files will 
never be required by HiggsBounds simultaneously. 
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file name data format 



BR_H_OP . dat * 


k, 


BRJij ss , 






BR Ji j c c , BRJi j bb , BRJi j t aut au , 






BR Ji j WW , BRJi j ZZ , BRJi j Zga, 






BRJijgaga.BRJijgg 


BR_H_NP . dat * 


k, 


BRJij invisible, some elements of BRJijhihi 






(row by row, without diagonal 






- see example) 


BR_t . dat * 


k, 


BR.tWpb, BR.tHpb 


BR_Hplus . dat * 


k, 


BR_Hpcs, BR_Hpcb, BR_Hptaunu 


CP_values.dat* 


k, 


CP.value 


additional . dat (optional) 


k, 





Table 13 

Names and data format of all HiggsBounds input Eles (part II). The right column shows 
the order of the input data arrays within one line of the input file. For the order within 
the arrays, see example, k is the line number. Note that the arrays CS_tev_pp_vbf _ratio, 
CS_tev_pp_tthj_ratio each appear in two different input files. However, these files will 
never be required by HiggsBounds simultaneously. 
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However, not all of the elements of the two dimensional arrays are required. Only the 
lower left triangle (including the diagonal) is required from the arrays g2hjhiZ and 
lepCS_hjhi_ratio, since they are symmetric, e.g. 

^g2hjhiZ(l, 1) g2hjhiZ(l.2) g2hjhiZ(l,3)^ 
g2hjhiZ(2, 1) g2hjhiZ(2,2) g2hjhiZ(2,3) 
^g2hjhiZ(3, 1) g2hjhiZ(3,2) g2hjhiZ(3,3)^ 

i.e. the elements in the input file should be written in the order 

g2hjhiZ(l,l), g2hjhiZ(2,l), g2hjhiZ(2 , 2) , g2hjhiZ(3, 1), 
g2hjhiZ(3,2), g2hjhiZ(3,3) . 

For the array BRJijhihi, only the off-diagonal components are required 

^BRJijhihi(l, 1) BRJijhihi(l, 2) BR_hjhihi(l, 3)^ 
BR_hjhihi(2, 1) BR_hjhihi(2, 2) BR_hjhihi(2, 3) 
^BRJijhihi(3, 1) BR_hjhihi(3, 2) BR_hjhihi(3, 3)^ 

since the diagonal elements are not physical quantities. Therefore, the elements should 
be written in the order 

BRJijhihi (1,2), BR_hjhihi(l,3), BR_hjhihi(2, 1), BR_hjhihi(2,3), 
BRJijhihi (3,1), BR_hjhihi(3,2) 

in the input file. 

The file additional.dat is optional. If it is included, it can have any number of 
columns greater than 1 (as for the previous files, the first entry on each line should 
be the line number). It is envisaged that this input file will be particularly use- 
ful when parameter scans are performed over a variable which is not required by 
HiggsBounds but helpful when plotting the results. For example, in the case of the 
MSSM, additional.dat could be used to store the values of tan/3. 

As in the subroutine version, the command line version of HiggsBounds expects 
a subset of the total list of input arrays, which depends on the chosen setting of 
whichinput. The maximal list of files used for each value of whichinput is given in 
Tab. 2. 

As discussed for the subroutine version, some of the arrays will not be relevant for 
some of the choices for whichanalyses. The command line version of HiggsBounds 
will consider the list of input files appropriate to the setting whichinput and then only 
attempt to read any of these input files if the value chosen for whichanalyses means 
that at least one of the arrays it contains will be directly used. Tab. 15 contains 
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a list of which input files are actually relevant to each value of whichanalyses. 
For example, if whichinput = 'hadr' , whichanalyses = ' LandT \ nHzero > and 
nHplus > then HiggsBounds requires the input files: 

MH_GainmaTot . dat, MHplus_GainmaTot . dat, CP_values . dat, 
BR_H_NP.dat, BR_H_OP.dat, BR_t.dat, BR_Hplus.dat, 

LEP_HZ_CS_ratios . dat, LEP_H_f f _CS_ratios . dat, LEP_2H_CS_ratios . dat, 
LEP_HpHm_CS_ratios . dat, TEV_lH_hadCS_ratios . dat . 

However, if whichinput = 'hadr', whichanalyses = 'onlyL', nHzero > and 
nHplus > HiggsBounds requires the input files: 

MH_GainmaTot.dat, MHplus_GammaTot.dat, CP_values.dat, 
BR_H_NP.dat, BR_H_OP.dat, BR_Hplus.dat, 

LEP_HZ_CS_ratios . dat, LEP_H_f f _CS_ratios . dat, LEP_2H_CS_ratios . dat, 
LEP_HpHm_CS_ratios .dat. 

As a third example, if whichinput = 'hadr', whichanalyses = 'onlyT', nHzero 
> and nHplus > HiggsBounds requires the input files: 

MH_GainmaTot.dat, MHplus_GainmaTot.dat, CP_values.dat, 
BR_H_NP.dat, BR_H_OP.dat, BR_t.dat, BR_Hplus.dat, 
TEV_lH_hadCS_ratios . dat . 

In each of these three examples, HiggsBounds will also read the file 

additional .dat 

if it exists. 

As for the subroutine version, if the user does not require processes involving a 
particular branching ratio or cross section ratio to be checked by HiggsBounds, that 
particular array can be filled with zeros. 



Output file format 

When the command line version of HiggsBounds is used, the output is written to the 
file <pref ix>HiggsBounds_results .dat. A sample of the output is shown in Fig. 1. 
The key to the process numbering is written to <pref ix>Key . dat. 

Examples 

The HiggsBounds package includes a full set of sample input files for the case uh = 
3, nH+ = 1, contained in the folder exainple_data. Each filename is prefixed with 
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whichinput — 'part', 


naar , 


) rt-F -FP ) 


Mn_bajnmal ot . aat 


Mn.bainmai ot . dat 


Mn.uammai ot . dat 


Mnplu.s_Gamniai ot .aat 


Mrlplus.bammai ot . dat 


Mrlplus.bammai ot . aat 


RR H MP Hal-* 


RR H MP Hal-* 


_^^p ^-4- (*) 

ex X . aax ^ ' 


r)n,_n_Ur . aat 


RD U nP /^-^-l- * 

DK_n_Ur . aat 


DD U "MP A'^+- * 

DK_n_iMr . aat 


riH_t . aat 


"D F) -4- J _ J- 

BH-t . aat 


"DD +- J _ -1- * 

DK_t . aat 


Dri_npiU.S . dat 


DK_npius . aat 


DK_npius . aat 


T TTP 117 pC T' o I- -i ri o 

j_iEjr^_tiZj_UD_i a.x.ios . aax. 


L.rjn_tiZj_>^o_i ax. XOS . aax. 


x<rjr^_npniii_uo_i axxos . aax 


jjE.Jr_ri_ii_L'0_raxios . aa.x 


T TTP U -F-F PQ -h -i o 

L.rjJr_n_II_L'0_rat lOS . aat 


additional . dat 


T 1713 OU /^C! -v '~> H — i ^ +- 


T T7P OU PC! A ^tn ^ +- 

LEir .zri-Oo-Tat los . aat 




J_iEjr^_tipniIl_L«o_i a.X.10S . QaX. 


LiEjr _npriiii_vvD_i ax.ios . aax. 




i Ej V _n_vj J eT._pai x.^-«o_x cLT. xos . ciaT. 


1 W ViaHfC; -i"a-|--ino Hal- (*) 

i Ej V _xn_QaaL/o_L ax XOS . aax ^ ' 




ihy Ji_ijex_parxuo_raT;ios . aat 


L.r_vaiues . aat 




ir.v_riw_parLUb_rat;ios . aat; 


additional . dat 




in.v_nZ_partOb_ratios . aat 






iE,v_n_VDi_n.aaCb_ratios . aat 






TEV_H_ttJiadCS_ratios . dat * 






CP_values.dat* 






additional . dat 







Table 14 

The list of possible input Eles for each value of whichinput. Note that some input files may 
not be relevant, depending on the values of whichanalyses, nHzero and nHplus. In this 
case, they are not required. See Tab. 15 and Tab. 16 for more details. 



HB_randomtest50points_. To run the command-line version of HiggsBounds with 
these files as input, use, for example, 

. /configure 
make 

./HiggsBounds LandT effC 3 1 ' example_data/HB_randomtest50points_ ' 

where the values of whichanalyses and whichinput can be varied as desired. The 
setting nHplus = can be used if the user does not wish to test the charged Higgs 
sector. E.g. 

./HiggsBounds LandT effC 3 ' example_data/HB_randomtest50points_ ' 
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name of input file values of whichanalyses 

which this file is relevant to 





LandT 


onlyL 


onlyT 


singH 


onlyP * 


jyin_bcLmmal ot . dat 


y 


y 


y 


y 


y 


rlnp±US_UaJllIIla. i OX. . QcLX. 


y 


y 


y 


y 


y 


ei I u . Qax. ^ ' 


y 


y 


y 


y 


y 


T T7D UV r\ci Ar\-^ 

J_lE.J^_^l^_L»o_J^axlos . aax 


y 


y 




y 


y 


7 170 U -F-F PQ -3 -1- -i o y-l-a-l-* 

JjEjJr_rl_ll_L'O_raXl0S . QaX 


y 


y 




y 


y 


T "FTP 0"H" PQ T' 1- -i /-i o 


y 


y 






y 


Lhir_ripnm_Oo_ratios . aat 


y 


y 




y 


y 


ihV_H_Ojet_partCb_ratios . aat 


y 




y 


y 


y 


iriV_h_ljet_partOb_ratios . aat 


y 




y 


y 


y 


ir.v_riw_pan;ob_raxios . aax 


y 




y 


y 


y 


iEiV_riZi_parLOo_rat;ios . aati 


y 




y 


y 


y 


iijv_n_VDi_LiaaL'0_raxios . ciax 


y 




y 


y 


y 


iiLv_n_xx_QacLL'0_i axios . aax 


y 




y 


y 


y 


irjV_in_LiaaL.o_raxios . aax ^ ^ 


y 




y 


y 


y 


BR_H_DP . dat 


y 


y 


y 


y 


y 


BR_H_NP . dat * 


y 


y 


y 


y 


y 


BR_t . dat * 


y 




y 


y 


y 


BR_Hplus . dat * 


y 


y 


y 


y 


y 


CP_values.dat* 


y 


y 


y 


y 


y 


additional.dat (optional) 


y 


y 


y 


y 


y 



Table 15 

hist of input Hies, specifying which values of whichanalyses each input file is relevant to 
(marked by 'y'). Note that some input files may not be relevant, depending on the values 
of whichinput, nHzero and nHplus. In this case, they are not required. See Tab. 14 and 
Tab. 16 for more details. 
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iicLiiit; oi input nit! 


UlUy ItJicVcilll 11 


uiiiy itJicVdiiL 11 




nHzero > 


nripius u 


nn_bajnmal ot . dat 


y 




Nnpius_v_raiiiniai or . aar 




y 




y 




T "FTP W PC; T'-a-!--; r.o 


y 




T FP H ff rc; -ra-|--inQ Hal-* 
J_il_jr^_n L J VjO L d. L ±Uo . Qd. L 


y 




T "FTP 0"H" PQ T' 1- -i /-i o 


y 




j_iEjr^_npniii_UD_i ax.ios . Qax. 




y 


i Hi V _n_uj eT._pai x.v-«D_x ax. xos . u.aT. 


y 




irjV_ri_ijet;_parrob_raxios . aax 


y 




iEiV_nW_parrL.o_raxios . aax 


y 




iEiV jiZ_parT;UD_raxios . aax 


y 




iiLV_ri_VDi_Liaaob_raxios . aax 


y 




iEiV_n_xx_QaaUo_raxios . aax 


y 




irjV_in_LiaaL.o_raxios . aax ^ ^ 


y 




RD TJ np rla-l- 

Drt_ti_ur^ . aax 


y 




OO TJ HTD J-N-l- * 

DHJUMr . aax 


y 




BH-t . dat 




y 


BR_Hplus . dat * 




y 


CP_values . dat * 


y 




additional . dat 


(optional) 


(optional) 



Table 16 



List of files, showing whicli relate to the neutral Higgs searches and which relate to the 
charged Higgs searches. Note that some input files may not be relevant, depending on the 
values of whichinput and whichanalyses. In this case, they are not required. See Tab. 14 
and Tab. 15 for more details. 
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# generated with HiggsBounds on 08.02.2011 at 11:18 

# settings: LandT, effC 
# 

# 
# 
# 
# 
# 
# 
# 
# 
# 
# 

# channel numbers used in this file 



column abbreviations 

line id of input 
masses of neutral Higgs boson 

masses of singly, positively charged Higgs boson masses 
scenario allowed flag (1: allowed, 0: excluded, -1: unphysical) 
most sensitive channel (see below) . chan=0 if no channel applies 
ratio [sig x BR] _model/ [sig x BR]_limit (<1: allowed, >1: excluded) 
number of Higgs bosons combined in most sensitive channel 
optional additional data stored in <pref ix>additional . dat (e.g. tan beta) 



Mh(i) 
Mhplus ( i ) 
HBresult 
Chan 

obsratio 

ncomb 

additional 



(ee)->(h3)Z->(b b)Z (LEP table 14b) 
(ee)->(hl)Z->(tau tau)Z (LEP table 14c) 
(pp)->W(hl)->l nu (b b) (CDF Note 9463) 
(pp)->h2->tau tau (arXiv: 0805 . 2491) 

(pp)->hl+... where hi is SM-like (arXiv: 0804. 3423 [hep-ex]) 



# 3 

# 4 

# 124 

# 134 

# 157 

# (for full list of processes, see Key. dat) 
# 

#cols 

# 



n 


Mh(l) 


Mh(2) 


Mh(3) 


Mhplus (1) 


HBresult 


Chan 


obsratio 




ncomb 


1 


359.121 


271.963 


134.929 


100.000 


1 


134 


0.212206E- 


-03 


1 


2 


75.0123 


92.8677 


71.9716 


100.000 


1 


4 


0.306172E- 


-01 


1 


3 


136.293 


345.483 


330.026 


100.000 


1 


124 


0.640713E- 


-01 


1 


4 


111.377 


220.765 


51.7469 


100.000 


1 


3 


0.162811 




1 


5 


186.131 


355.002 


146.448 


100.000 





157 


15.2354 




1 



additional (1) 

0.246862 
0.714964 
0.434594 
0.727173 
0.230522 



Fig. 1. Sample output file (written to <pref ix>HiggsBoundsJlesults . datj 

5.5 Online version 



The online version can be reached via the website ^ 

http : / /pro j ects . hepf orge . org/higgsbounds/ . 

It allows the user to select the required number of neutral Higgs and charged Higgs 
bosons and then generates a html form accordingly. The values of whichinput and 
whichanalyses can be chosen and the appropriate theoretical input entered. Higgs- 
Bounds will then be called with these settings and the result outputted to screen. 
The online version contains the additional feature that it notifies the user about the 
processes with the second and third highest statistical sensitivities and the values 
of obsratio for these processes. This is designed to give guidance to the user who, 
for example, wishes to find an excluded region iteratively by adjusting the input 
quantities. 

The website also contains a selection of pre-filled html forms as examples, including 
entries for the SM, the fermiophobic Higgs model and the MSSM with real and 
complex parameters. 

* There is also the possibility of copying and pasting the results from the online 
version of FeynHiggs (the 'FeynHiggs User Control Center' [111]) into a box on the 



^ The former website www . ippp . dur . ac . ulc/HiggsBounds should redirect to the new one. 
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HiggsBounds website for the CP-conserving MSSM. This text is then converted to 
a filled-in HiggsBounds input html form (which can be further edited by the user if 
necessary) that can then be immediately submitted to HiggsBounds. 



6 Physics Applications 



6.1 Model scenarios with invisible Higgs decay modes 



The decay of a Higgs boson into stable, weakly interacting neutral particles (invisible 
Higgs decay) is allowed in a wide variety of models. In the R-parity conserving MSSM 
with a neutralino Xi the lightest supersymmetric particle (LSP), the decay — )■ 
XiXi is dominant in some scenarios [112,113]. Invisible Higgs decays can also be 
important in models with hidden sectors that couple to the Higgs sector [114], models 
where a Higgs boson can decay into a pair of Majorons [115-117], models with 4th 
generation neutrinos [118], non-linear supersymmetric models in which the Higgs 
boson can decay into a Goldstino and a neutralino [119] and some extra dimension 
models [120-125]. 

We consider here a simple model in which we vary the branching ratio of the invisible 
Higgs decay mode, while keeping all other cross sections and decay widths fixed to 
SM values and assuming that the narrow width approximation holds. This is achieved 
with the following program: 

j3|c*;|cHc**3|(Hc************* ************************* 

program example.invisible 
implicit none 



integer 
integer 
integer 



nH , nHplus 

HBresult , chan , ncombined 
p,q,pend,qend 

obsratio 
mass 

SMGammaTotal 



double precision 
double precision 
double precision 
! HB input: 

double precision :: Mh,GammaTotal_hj , & 
& g2hj ss_s ,g2hjss_p,g2hj cc_s , g2hj cc_p, & 

& g2hj bb_s , g2hjbb_p , g2hj t opt op_s , g2h j toptop_p , & 

& g2hjtautau_s,g2hjtautau_p, & 

& g2hjWW,g2hjZZ,g2hjZga, k 

& g2hjgaga,g2hjgg,g2hjhiZ_nHbynH, & 

& BR_hj invisible, BR_hjhihi_nHbynH 

! HB functions: 

double precision : : SMGamma_h , SMBR_Hbb 

nH=l ! only one neutral Higgs in model 

nHplus=0 ! no charged Higgs in model 

call initialize_HiggsBounds(nH, nHplus, 'LandT') 

pend=201 

qend=201 

open (10, file=' example-invisible-result s .dat ' ) 



HB input 
HB output 
used in do loops 
HB output 
Higgs mass 

SM total decay width for a SM Higgs at this mass 
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do p=l,pend 

Mh= 50.0D0 + dble(p-l)*200.0D0/dble(pend-l) 
do q=l,qend 
BR_hj invisible= dble (q-1) *1 . ODO/dble (qend-1) 

SMGanraiaTotal = SMGamma_h(Mh) 

if (SMGammaTotal.lt. O.ODO)stop'negative result from SMGamma_h(Mh) ' 

if (abs ( 1 . ODO-BR_hj invisible) . gt . . ODD) then 

GammaTotal_hj = SMGammaTotal/ (1 . ODO-BR_hjinvisible) 
else 

GaiiimaTotal_hj = 1.0D8 ! irrelevant if higgs->invisible is the only decay mode 
end if 



set normalised couplings to SM values i.e. 
all except pseudoscalar couplings set to 1 
pseudoscalar couplings set to 



g2hjss_s 


= 1 


ODD; 


g2hjss_p 


= 


ODO 


g2hjcc_s 


= 1 


ODD; 


g2hjcc_p 


= 


ODO 


g2hjbb_s 


= 1 


ODD; 


g2hjbb_p 


= 


ODO 


g2hjtoptop_s 


= 1 


ODO; 


g2hjtoptop_p 


= 


ODO 


g2hjtautau_s 


= 1 


ODO; 


g2hjtautau_p 


= 


ODO 



g2hjWW =1.0D0 
g2hjZZ =1.0D0 
g2hjZga =1.0D0 
g2hjgaga=1.0D0 
g2hjgg =1.0D0 
g2hj hiZ_nHbynH=0 . ODO 



! only one Higgs in model, so no Higgs cascade decay branching ratio 
BR_hjhihi_nHbynH=0 . ODO 

call HiggsBounds_neutral_input_ef f C( Mh,GammaTotal_hj , & 



& g2hj ss_s ,g2hjss_p,g2hj cc_s ,g2hj cc_p, & 

& g2hjbb_s ,g2hjbb_p,g2hjtoptop_s ,g2hjtoptop_p,& 

& g2hjtautau_s ,g2hjtautau_p , k 

& g2hjWW,g2hjZZ,g2hjZga, k 

k g2hjgaga,g2hjgg,g2hjhiZ_nHbynH, & 

k BR_hj invisible, BR_hjhihi_nHbynH ) 



call run_HiggsBounds ( HBresult , chan, obsratio, ncombined ) 

write (10 , *) Mh,BR_hj invisible , HBresult , chan, obsratio 

enddo 
enddo 

call f inish_HiggsBounds 
close(lO) 

end program example_invisible 

which can be compiled using, for example, 

gfortran . /example_invisible . f 90 -o . /example_invisible . exe -L<HBdirectorypath> -IHB 



6.1.1 Exclusion 



We show the results from the code displayed above in Fig. 2. In the left panel, we can 
see that all values of the invisible Higgs branching ratio are excluded almost up to 
the kinematical limit of LEP. In the right panel, we can see that LEP Higgsstrahlung 
processes have the highest statistical sensitivity in this region, and therefore were 
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exclusions at 95% CL analysis with highest stat. sens. 




50 100 150 200 250 50 100 150 200 250 



Mn/GeV Mn/GeV 



Fig. 2. Parameter space for our toy model, in which we vary the branching ratio of the 
Higgs decay into invisible particles and the Higgs mass. Higgs production cross sections 
and all other Higgs decay widths take SM values. Left: parameter points excluded at 95% 
CL. by Higgs searches at LEP and the Tevatron. Green (dark grey) = excluded, White = 
unexcluded. 

Right: Channel with the highest statistical sensitivity at each parameter point. 

U = e-e+ HZ bbZ, Ref [2] (LHWG) 

U = e'e^ HZ (invisible) Z, Ref [31] (OPAL) 

U = e-e+ HZ (invisible) Z, Ref [30] (L3) 

= e-e+ HZ (invisible) Z, Ref [28] (LHWG) 
I =pp^H ^W+W-, Ref [5] (TEVNPHWG) 
U=pp^ H + ...^ ... where H is SM-like, Ref [85] (TEVNPHWG) 



used to obtain this exclusion. For lower values of invisible Higgs branching ratio, 
the decay mode H ^ bb was used, whereas at higher values of BR(if — T-invisible), 
Higgsstrahlung topologies which explicitly involve the H — j-invisible decay mode (i.e. 
searches for a Z-boson plus large missing momentum) were used. The LEP Higgs 
Working Group combined invisible Higgs analysis [28] does not cover the region Mh < 
90 GeV, and therefore the OPAL [31] and L3 [30] analyses are also required. 

We can also apply the results from the Tevatron SM Higgs searches directly to this 
model, since all the relevant cross sections and branching ratios differ from their SM 
equivalents by a common factor (i.e. the parameter points all pass the SM-likeness 
test for the Tevatron SM Higgs analyses, as described in Sect. 4). Therefore, the 
exclusion at 95% CL. of a SM Higgs boson between 158 GeV and 175 GeV provided 
by Ref. [85] translates to a wedge-shaped excluded region in the parameter space of 
our toy model, extending up to BR(H — j'invisible)=0.32. 
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Fig. 3. Subsections of the parameter space for our toy model, in which we vary the branching 
ratio of Higgs into invisible particles and the Higgs mass, obsratio is the ratio of the theo- 
retical cross section to the observed cross section for the channel with the highest statistical 
sensitivity, as calculated by the program example_invisible, plotted when obsratio > 1 
(left) or when obsratio < 1 (right). Left: we interpret i/obsratio as the maximum value 
of 

9hzz which is consistent with the parameter point being excluded at 95% C.L. Right: 
we interpret 1/obsratio as the factor /exci by which the measured SM cross section limit 
must be reduced in order to exclude this parameter point at 95% C.L. (in this plot, red 
denotes /cxci > W). 

6.1.2 Interpreting the variable obsratio 

We can also obtain interesting information about our model from the variable obsratio 
outputted by the program example_invisible. Recall that obsratio is the ratio of 
the theoretical cross section to the observed cross section for the channel with the 
highest statistical sensitivity at each parameter point Qmodei(^o)/Qobs(^o), as dis- 
cussed in Sect. 5. Roughly speaking, this variable gives an idea of 'how strong' the 
exclusion of a particular parameter point is. We shall now demonstrate two possible 
ways in which this variable can be used to obtain additional insights into our toy 
model. We will first discuss the excluded low mass region (with obsratio > 1) and 
then we will discuss the unexcluded intermediate mass region (with obsratio < 1). 

In models with extended Higgs sectors, such as the MSSM, the Higgs couplings to 
Z-bosons frequently obey the sum rule Y.i=i Qh zz — 1; therefore the hi-Z-Z cou- 
plings in these models are suppressed with respect to the SM coupling. It is therefore 
interesting to consider the extent to which the parameter space of our toy model can 
still be excluded if we suppress the H-Z-Z coupling. We can interpret the inverse of 
the obsratio calculated by the program example_invisible as the maximum gl^zz 
which will still allow the parameter point to be excluded at 95% by the LEP Higgs 
search results (considering Higgsstrahlung processes only). 1/obsratio is plotted in 
Fig. 3(left) in the low Higgs mass region, for obsratio > 1. In this plot, we can 
see that, for example, a Higgs boson with Mh = 95 GeV and BR(if ^■invisible)=0.5 
(and all other decay widths as in SM) can not be excluded by the LEP Higgsstrahlung 
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searches if gjjzz < 



We shall now look at 1/obsratio in the unexcluded part of the intermediate Higgs 
mass range, where the Tevatron SM Higgs analysis Ref. [85] has the highest statistical 
sensitivity. This is shown in Fig. 3(right), for obsratio < 1. In this mass range, it is 
interesting to interpret 1/obsratio as the factor /exci by which the current SM cross 
section limit must be lowered in order to obtain a 95% exclusion in our toy model. 
Therefore, if, for example, a future Tevatron analysis can lower the current limit by 
10% at Mh = 165 GeV, this will allow BR{H -^invisible) up to 0.38 to be excluded 
in our toy model at this Higgs mass. In this way, we can get an impression of which 
areas of parameter space are likely to be testable in the near future and which areas 
will require much more data. 

6.2 Constraints on the Randall- Sundrum scalar sector 

The Randall- Sundrum (RS) Model, RSI, considers spacetime a slice of 5d anti-de- 
Sitter space with two boundaries, the IR brane (our 4d spacetime) and the UV brane 
[126]. This configuration yields an explanation for the hierarchy problem, i.e. why 
the mass scale of electroweak physics is so much smaller than the Planck scale Mpi. 
The spacetime metric can be written as 

ds^ = e-^^'^yr^^^dx^'dx" - rldy^ , y G [0, tt] , 

with (yU G 0, 1,2,3) coordinates of the usual 4d spacetime and y the coordinate 
of the fourth spatial direction. The quantities k, are 0{Mpi) with kr^ ~ 12. The 
"little hierarchy" between k and can be generated and stabilized naturally [127]. 
The hierarchy between electroweak and Planck scale is a consequence of the "warped" 
metric: mass parameters in the fundamental 5d model mo appear in our visible space 
as 

m = moe"^'""'' ^ molO"^*^ . 

Hence, all mass scales of the model can be assumed to be of the order of the Planck 
scale. Many modifications to the original RS Model, concerning which fields should be 
localised on the branes, have been considered necessary in order to meet electroweak 
precision and flavour constraints [128]. However, the Higgs field needs to be localised 
on or near the IR brane in order to explain the hierarchy problem. Therefore, the RS 
scalar sector is a rather robust prediction of RS models. 

6.2.1 The Higgs-radion system 

Fluctuations of the compactification "radius" Vc correspond to a scalar degree of 
freedom in addition to the 4d metric fluctuations. As a consequence of stabilising rc, 
this scalar acquires a vacuum expectation value (VEV) and appears as a massive 
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scalar (the radion (f)oix)) in the spectrum [127,129]. Higgs-radion mixing may occur 
via the curvature-scalar interaction [130] 



^ = -^V-C/ind R{9md) , 

where gmd{4>o{^)y ' ' ' ) is the induced 4d metric on the IR brane, R the Ricci scalar, 
and $ the Higgs doublet. Hence, in general, what would be for ^ = the physical 
Higgs degree of freedom ho and the radion 0o Hiix to form two mass eigenstates h 
and (p. One arrives at the fields Hq and (po after rescaling to canonical normalisation 
on the IR brane. The part of the Lagrangian bilinear in ho and (po reads [131,132]: 

= -^(1 + 67'O0on0o - IhoDho - 67^000/^0 - ^mj„0o - ^^L^o • (37) 

Diagonalisation and canonical normalisation of the bilinear terms can be obtained 
for the physical fields h and via the identification: 



a(f) + bh , ho = dh + ccj) , (38) 



with 



a = — , b = ^" ^ , c = sin 6* + cos 9 , d = cos 9 — sin 9 , (39) 
Zj Zj Zj z 

and 

t^^29 = — 2 ^2 \a,2 2W ^' = 1 + 6^7' 1-60, 7=-r-. 40 

In this formulae, v (~ 246 GeV) is the Higgs VEV. The radion VEV A,^ is a free pa- 
rameter of the model. The masses and couplings to known matter and gauge bosons 
of the Higgs-radion sector are determined by choosing values for the four parame- 
ters A^, the curvature-Higgs coupling parameter ^, and the masses nih and of 
the physical fields. However, those parameter choices are subject to theoretical con- 
straints, which come from requiring positive kinetic energy and mass terms in the 
Lagrangian [131,132]: 



> , max 



2 ' 2 



> 1 + - (VT^ +l-x) , (41) 



with X = (l-|-36.^^7^/Z^) ^. The Lagrangian parameters m^o and m^o are determined 
for a given valid parameter scenario by 

X Z"^ 
mlo = -(S + [-]A) , m^o = — (S - [+]A) , if mt > m^ [mh < m^] , 

I 4 

with S = ml + m^ , A = W 'mfjn?^ ■ 
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In the following, we do not consider scenarios where the physical scalars can decay 
into Kaluza-Klein (KK) excitations of the graviton, i.e. implicitly we assume that 
the mass mi of this excitation is large enough to forbid such decay channels. This 
puts a lower bound on the curvature parameter k. 

Providing values for A^, ^, nih and for a valid scenario, the quantities a, b, c, and 
d defined in Eq. (38) are determined. This, in turn, determines the effective couplings 
of the two physical scalars h and to gauge bosons and matter normalised to the 
couplings of a SM Higgs boson with the same mass [120,132]: 



nh ^ //) 
nHsM ^ //) 



T{HsM gg) 
r(0 ^ gg) 
r(iJsM gg) 

^{h 77) 
T{HsM 77) 

r(0 ^ 77) 



r{HsM w+w- 

T{(f) W+W-) 



2 1266, 



T{h ZZ) 
T{HsM ^ ZZ) 

T{(j)^ ZZ) 
T{HsM ^ ZZ) 



V 



w+w-\ 

ib+^d)Fy2iTt)\' 



A2 



\Fy2iT!^)\' 
2ab,-{a + ^c)Fy,{Tt)\' 



2 \3b{b2 



\Fi/2{rr)\' 



d + 








V 


c + 


it"- 





V 



{b+^d){AFy,irl^) + 3F,iTm' 



\3a{b2 + by) 



|4Fi/2(ri'^) + 3Fi(r^) 



a + ^c)(4Fv2(rn + 3Fi(r^))P 



T{HsM^ll) Al |4Fi/2(r/) + 3Fi(r^ 

with 62 = 19/6, 63 = 7, 6y = -41/6, Ti'^'^l = 4m?/m2[^,, r^^l = Aml,/ml^, 



and 



FiMr) = -2r(l + (1 - r)/(r)) , F,{t) = 2 + 3r + 3r(2 - r)/(r) , 



/( 



arcsin^fr ^^'^) 



T 



In 



i+\/T^ 
l-\/T^ 



T > 1 
T > 1 



Like for the Higgs, the radion couplings to massive fermions and gauge bosons are 
proportional to mass, but e.g. the couplings 0o bb and 0o 77 are suppressed while 
00 gg is enhanced with respect to the SM Higgs boson. In order to consider all possible 
search channels for the Randall- Sundrum model with HiggsBounds, knowledge of the 
branching ratios for the decays /i — )■ 00 and — )■ /i/i is needed as well . However, 
those channels turn out not to be the most significant ones in the examples we study 
below. 
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Formulas for those branching ratios can be obtained from [132]. 
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6.2.2 Constraints 



We chose here a rather extreme scenario with as low as 1 TeV, the phenomenology 
of which has been pioneered in [132]. This scenario may already have been excluded 
by other observations. In particular, our assumption that mi is high enough such 
that decays oi h oi cf) into graviton KK excitations are impossible, entails a rather 
high curvature scale k in conflict e.g. with direct searches for graviton KK excitations 
in the Randall-Sundrum model [133]. However, it serves demonstration purpose well. 
To our knowledge, this is the first time that Tevatron Higgs search results have been 
recast in order to provide constraints on the scalar sector of the Randall-Sundrum 
model. A more thorough phenomenological study will appear elsewhere [134]. 

Fig. 4 shows for the RS model the excluded region (left panel) and the search channel 
with highest sensitivity (right panel) in the — plane, while the two other free 
parameters in this model are set to A,^ = 1 TeV and ^ = 1/6. Note the region of 
parameter space inaccessible due to the theoretical constraints in Eq. (41) shown in 
Fig. 4 (left panel). This scenario shows slight Higgs-radion mixing, i.e. h behaves 
mainly like the SM Higgs and mainly like the unmixed radion. This is reflected 
by the LEP exclusion of values up to almost the SM mass limit via the process 
e'^e~ — hZ, h hh (see Fig. 4, right panel) almost independent of m^. For larger 
values of m/j, there is a LEP excluded region via e^e~ — > </)Z, </> — hadrons for 
below about 50 GeV. In this region, the enhanced gg coupling and the suppressed hb 
coupling of the radion-like renders this channel the most sensitive one. The exclusion 
of parameter regions by Tevatron searches is mainly due to the model independent 
cross section limit on the process pp — ?■ single S, 5* — )■ W^W~ — )■ Ivlv [5] for a 
scalar S. Indeed, for scalar masses mh and m<^ above around 120 GeV in the displayed 
parameter region, this analysis has the highest statistical sensitivity (see Fig. 4, right 
panel). Furthermore, the highest exclusion power of this analysis is for scalar masses 
around ^ 165 GeV. The right panel of Fig. 4 reveals that sizable regions where 
either or have values close to 165 GeV are excluded by this analysis with S = h 
or S = (j). For instance for nih = 300 GeV, the excluded range of values due to 
non-observation of the radion-like scalar is wider than for a SM-like Higgs because of 
the net enhancement of the signal cross section a{pp gg ^ (p) x BR(</) — )■ iy"'"iy~) 
with respect to the SM. 

Fig. 5 shows the excluded Randall-Sundrum parameter space also for A^ = 1 TeV 
but as function of ^ and with fixed nih = 120 GeV. The parameter exclusion for 
above 120 GeV is entirely due to the application of the Tevatron analysis [5] to the 
radion-like 0, i.e. due to the search channel pp — )■ single (f),(j) ^ W~^W~ — t- lulu (see 
Fig. 5, right panel). In this region, the effective gg(f) coupling is rising with falling ^, 
while the effective WWcf) coupling goes through a minimum located near ,^ ^ 0. This 
explains, on the one hand, the general trend of an increasing excluded m^^-interval 
with falling ^ values and, on the other hand, the unexcluded funnel region near ^ k, {). 



Fig. 5 can be seen as an update of Fig. 18 in [132]. 
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Fig. 4. Excluded region in the rrih — plane of the scalar sector of the original Randall- 
Sundrum model (RSI) for = 1 TeV and ^ = 1/6 (left panel) and most sensitive search 
channel (right panel) using' HiggsBounds 2.0.0. 





Randall-Sundrum Model (RSI) 
A,j = 1 TeV, mk = 120 GeV 

most sensitive analysis: 

■ : e+e- ^ hZJi 66 
■: e+e- (jiZ, bb 
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■ : pp -» single h,h^ W*W~ 

■ : pp — f single if),(j) ^ W^W~ 



Fig. 5. Excluded region in the ^ — plane of the scalar sector of the original Randall-Sun- 
drum model (RSI) for = 1 TeV and = 120 GeV (left panel, see Fig. 4 (left panel) 
for the key) and most sensitive search channel (right panel) using HiggsBounds 2.0.0. 

For below 120 GeV, another region is excluded by applying the Tevatron analysis 
[5] — this time to the Higgs-like h. This analysis application becomes the most 
sensitive one because in the lower left part of the allowed parameter space, both the 
effective ggh and WWh coupling rise sufficiently with falling ^. The LEP excluded 
region consists of one region where the most sensitive channel is e~^e— — (pZ, cf) — 
hadrons [35-38] (for roughly below 50 GeV in Fig. 5, right panel), and another 
one where it is e^e > (pZ, cf) -^hh [2]. 



7 Summary 



We have presented the code HiggsBounds 2.0.0 which allows to test neutral and 
charged Higgs sectors of arbitrary models against the current exclusion bounds from 
the Higgs searches at LEP and the Tevatron. The model predictions which are re- 
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quired as input to the program, and in particular the three variants of how this 
information can be provided by the user, are described in detail. Version 2.0.0 repre- 
sents a significant upgrade of the code since its first release. As a major extension, the 
code allows now the predictions for (singly) charged Higgs bosons to be confronted 
with LEP and Tevatron searches. Furthermore, the newly included analyses contain 
LEP searches for neutral Higgs bosons (H) decaying invisibly or into (non flavour 
tagged) hadrons as well as decay-mode independent searches, LEP searches via the 
production modes t^t~H and hbH, and Tevatron searches via tiH . Many updated 
Tevatron analyses have replaced previously included ones. In particular, all Tevatron 
results presented at the ICHEP'IO are included in HiggsBounds 2.0.0. We presented 
a full list of the analyses (with references) incorporated in the program and the im- 
plemented conditions for Higgs bosons to be met in order to be applicable to some of 
the analyses. We provide the complete operating instructions for HiggsBounds2.0.0 
with examples, indicating new or extended features with respect to version 1.0.0 of 
the code. 

We have furthermore presented two new physics applications of the program. One 
application considers a SM-like Higgs boson where an extra Higgs decay channel to 
invisible particles is open. We find, for instance, that the current Tevatron searches 
for a SM-like Higgs boson yield an exclusion interval around tuh ^ 2my/ even for 
the case of a Higgs boson which decays invisibly with up to 32% branching ratio. As 
a further application we have investigated the scalar sector of the Randall- Sundrum 
model. Our results show, to our knowledge, for the first time the impact of Tevatron 
Higgs search results on this model. For the parameter scenario we study here, we 
find that once the Tevatron results are taken into account in addition to the LEP 
constraints that had been considered before in [132], large portions of the parameter 
space are excluded as a consequence of the model-independent Tevatron Higgs search 
via the decay H W^W~ [5]. 

The functionality of HiggsBounds 2.0.0 is such that upcoming Higgs search limits 
from the LHC can easily be incorporated. The bounds from Higgs searches at LEP, 
the Tevatron and the LHC itself will provide an important source of information for 
testing model interpretations of possible Higgs-like signals, which will hopefully soon 
emerge from the LHC Higgs searches. 
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